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E. Jurcisinova 1,2 , and M. Jurcisin 1,3 
1 Institute of Experimental Physics, Slovak Academy of Sciences, Watsonova 47, 040 01 Kosice, Slovakia 
2 Laboratory of Information Technologies, Joint Institute for Nuclear Research, 141 980 Dubna, Moscow Region, Russia 

3 N.N. Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 141 980 Dubna, Moscow Region, Russia 

(Dated: February 8, 2008) 

The influence of uniaxial small-scale anisotropy on the stability of the scaling regimes and on the 
anomalous scaling of the structure functions of a passive scalar advected by a Gaussian solenoidal 
velocity field with finite correlation time is investigated by the field theoretic renormalization group 
and operator product expansion within one-loop approximation. Possible scaling regimes are found 
and classified in the plane of exponents e — n, where e characterizes the energy spectrum of the 
velocity field in the inertial range E oc fc 1_2e , and n is related to the correlation time of the velocity 
field at the wave number k which is scaled as k~ 2+v . It is shown that the presence of anisotropy 
does not disturb the stability of the infrared fixed points of the renormalization group equations 
which are directly related to the corresponding scaling regimes. The influence of anisotropy on the 
anomalous scaling of the structure functions of the passive scalar field is studied as a function of 
the fixed point value of the parameter u which represents the ratio of turnover time of scalar field 
and velocity correlation time. It is shown that the corresponding one-loop anomalous dimensions, 
which are the same (universal) for all particular models with concrete value of u in the isotropic 
case, are different (nonuniversal) in the case with the presence of small-scale anisotropy and they are 
continuous functions of the anisotropy parameters, as well as the parameter it. The dependence of 
the anomalous dimensions on the anisotropy parameters of two special limits of the general model, 
namely, the rapid-change model and the frozen velocity field model, are found when u — * 00 and 
u — > 0, respectively. 

PACS numbers: 47.27.-i, 47.10.+g, 05.10.Cc 
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I. INTRODUCTION 

One of the last unsolved problem in the framework of 
classical physics still remains the theoretical understand- 
ing of turbulence. Within one part of the comprehen- 
sive concept of turbulence, namely, fully developed tur- 
bulence, one of the most interesting and still open ques- 
tion is the theoretical explanation and understanding of 
the possible deviations from the classical phenomenolog- 
ical Kolmogorov-Obukhov theory which are suggested by 
both natural, as well as numerical experiments [J 0, La, 0] • 
Such a behavior is contained in concepts intermittency 
and anomalous scaling. During the last two decades 
this problem was intensively studied within the scope of 
the models of passively advected scalar field (concentra- 
tion of an admixture, or temperature are examples) by 
a velocity field with prescribed Gaussian statistics. The 
reason is twofold. First, it is well known that the de- 
viation from the classical Kolmogorov-Obukhov theory 
is even more strongly noticeable for passively advected 
scalar field then for the velocity field itself, see, e.g., 
Ref. 0, i, H, 0, 0, 1 E, Q II 0, and second, the 
problem of passive advection of a scalar or vector field 
is considerably easier for theoretical investigation than 
the original problem of anomalous scaling in the frame- 
work of Navier-Stokes velocity field. On the other hand, 
these simplified models reproduce many of the anomalous 
features of genuine turbulent heat or mass transport ob- 
served in experiments. Thus, on one hand, the theoretical 



study of the models of passive scalar (or also vector) ad- 
vection can be treated as the first step on the long way of 
the investigation of intermittency and anomalous scaling 
in fully developed turbulence but, on the other hand, the 
problem of advection has its own practical importance 
(see, e.g., Ref. and references cited therein). 

The central role in the studies of passive advection 
was played by a simple model of passive scalar quan- 
tity advected by a random Gaussian velocity field, white 
in time and self-similar in space, the so-called Kraichnan 
rapid-change model jl4j |. Namely, in the framework of 
the rapid-change model, for the first time, the anoma- 
lous scaling was established on the basis of a microscopic 
model [15| and corresponding anomalous exponents were 
calculated within controlled approximations [l6|, [l7| (see 
also survey paper [l8| and references cited therein). 

An effective method for investigation of self-similar 
scalin g be havior is the renormalization group (RG) tech- 
nique |ljj HH • It plays crucial role in the explanation 
of the origin of critical scaling in the theory of critical 
phenomena, as well as it allows to calculate some univer- 
sal quantities (e.g., critical dimensions). The RG tech- 
nique can be also used in the theory of fully developed 
turbulence and related problems [2l|, [2^, [H, H3] (passive 
advection is an example). It is important to note that 
there are many different RG methods, with the same idea 
but with technical differences, but perhaps the most for- 
malized one is the so-called "quantum field theory" RG 
which is also known as "field theoretic" RG [2l[. It is 
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based on the standard renormalization procedure, i.e., 
on the elimination of ultraviolet (UV) divergences. 

In Ref. [25| the field theoretic RG and operator-product 
expansion (OPE) was used in the systematic investiga- 
tion of the Kraichnan's rapid-change model, where it 
was shown that within the field theoretic RG approach 
the anomalous scaling is related to the existence in the 
model of the composite operators with negative critical 
dimensions in the OPE which are usually termed as dan- 
gerous operators (see, e.g., [HI, US 0] for details). In 
Ref. [25I ] the anomalous exponents were calculated to or- 
der e 2 (two-loop approximation) within the e expansion 
where parameter e describes a given equal-time pair cor- 
relation function of velocity field (see subsequent section) 
but quite early after this important work papers (26j have 
appeared where the power of the field theoretic RG was 
fully demonstrated, namely, the anomalous exponents of 
the Kraichnan model were calculated to order e 3 (three- 
loop approximation) within the e expansion. This result 
was not achieved by any other method yet and as far as 
we know this is the only known three-loop result in fully 
developed turbulence and related problems at all. 

Afterwards, various descendants of the Kraichnan 
model, namely, models with inclusion of small scale 
anisotropy [13], com pres sibi lity l28l . I29I ] . finite correlation 
time of velocity field [30l.l3lll32l.l33j]. and helicity [H[ were 
studied by field theoretic approach. Moreover, advection 
of passive vector field by Gaussian self-similar velocity 
field (with and without large and small scale anisotropy, 
pressure, compressibility, and finite correlation time) has 
been also investigated and all possible asymptotic scaling 
regimes and cross-over among them have been classified 
and anomalous scaling was investigated [35]. General 
conclusion of all these investigations is that the anoma- 
lous scaling, which is the most intriguing and important 
feature of the Kraichnan rapid change model, remains 
valid for all generalized models. 

The Kraichnan model works with white in time (5 cor- 
related in time) and self-similar in space Gaussian statis- 
tics of the velocity field. In Ref. [3(| the field theoretic 
RG technique and OPE method was applied in the anal- 
ysis of more general model of passively advected scalar 
field by a self-similar Gaussian velocity field with finite 
correlation time first proposed in Ref. This model 
contains the Kraichnan model as a special limit case (see 
next section). Maybe the most interesting conclusion 
from the view of anomalous scaling analysis obtained 
in Ref. [3(j is that within the one-loop approximation 
the anomalous behavior of all particular models of the 
general one (the Kraichnan model is an example) is the 
same, i.e., the corresponding critical dimensions associ- 
ated with needed composite operators within the OPE 
are the same. This conclusion is held in isotropic model, 
as well as in the model with large-scale anisotropy with 
incompressible (solenoidal) velocity field. This univer- 
sality of the anomalous behavior is destroyed, e.g., by 
the assumption that velocity field is nonsolenoidal as was 
shown in Ref. (3lj or by the assumption of the presence 



of small-scale anisotropy of the velocity field what will 
be demonstrated explicitly in present work. But first let 
us motivate the importance of such investigations. 

In Ref. [13] the field theoretic RG and OPE were ap- 
plied to the rapid change model of passive scalar advected 
by Gaussian strongly anisotropic velocity field where the 
anomalous exponents of the structure functions were cal- 
culated to the first order in e expansion. It was shown 
that in the presence of small-scale anisotropy the corre- 
sponding exponents are nonuniversal, i.e., they are func- 
tions of the anisotropy parameters, and they form the 
hierarchy with the leading exponent related to the most 
"isotropic" operator. The importance of these investi- 
gations is dictated by the question of the influence of 
anisotropy on inertial-range behavior of passively ad- 
vected fields [H M, m Er 

as the velocity field itself [4 



MS SO, 53 > as wel1 
44j (see also the survey 



paper [4f| and references cited therein, as well as recent 
astrophysical investigations, e.g, in Refs. [46l . H?))- On 
one hand, it was shown that for the even structure (or 
correlation) functions the exponents which describe the 
inertial-range scaling exhibit universality and they are 
ordered hierarchically in respect to degree of anisotropy 
with leading contribution given by the exponent from 
the isotropic shell but, on the other hand, the survival 
of the anisotropy in the inertial-range is demonstrated 
by the behavior of the odd structure functions, namely, 
the so-called skewness factor decreases down the scales 
slower than expected earlier in accordance with the clas- 
sical Kolmogorov-Obukhov theory. 

Let us describe briefly the solution of the problem in 
the framework of the field theoretic approach [U [23|, [24[ ■ 
It can be divided into two main stages. On the first stage 
the multiplicative renormalizability of the corresponding 
field theoretic model is demonstrated and the differential 
RG equations for its correlation functions are obtained. 
The asymptotic behavior of the latter on their ultraviolet 
argument (r/£) for r 3> I and any fixed (r/L) is given 
by infrared stable fixed points of those equations. Here I 
and L are inner (ultraviolet) and outer (infrared) scales 
(lengths). It involves some "scaling functions" of the in- 
frared argument (r/L), whose form is not determined by 
the RG equations. On the second stage, their behavior 
at r C I is found from the OPE within the framework 
of the general solution of the RG equations. There, the 
crucial role is played by the critical dimensions of var- 
ious composite operators, which give rise to an infinite 
family of independent aforementioned scaling exponents 
(and hence to multiscaling). 

In Ref. [10] the problem of a passive scalar advected 
by Gaussian self-similar velocity field with finite cor- 
relation time [48[ was studied by field theoretic RG 
method. There, the systematic study of the possible scal- 
ing regimes and anomalous behavior was present at one- 
loop level. The two-loop corrections to the anomalous 
exponents were obtained in Ref. [32] . It was shown that 
the anomalous exponents are nonuniversal as a result of 
their dependence on a dimensionless parameter u, the ra- 
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tio of the velocity correlation time, and turnover time of 
scalar field. 

In what follows we shall continue with the investiga- 
tion of this model from the point of view of the influence 
of the uniaxial small-scale anisotropy on the anomalous 
scaling of the single-time structure functions. In con- 
tradistinction with the studies of 30] , where the velocity 
was isotropic and the large-scale anisotropy was intro- 
duced by the imposed linear mean gradient, the uniaxial 
anisotropy in our model persists for all scales, leading 
to nonuniversality of the anomalous exponents through 
their dependence on the anisotropy parameters and ra- 
tio of characteristic time scales. It can be consider as 
an additional step to the construction of a more realistic 
model of anisotropic passive advection. 

The aim of the present paper is twofold. First of all 
we shall find the dependence of the anomalous expo- 
nents on the anisotropy parameters of the model and 
on the parameter u, therefore we shall be able to answer 
the question whether the system with finite time corre- 
lations of the velocity field with presence of small-scale 
anisotropy is more anomalous, i.e., whether the corre- 
sponding critical dimensions are less than those of the 
Kraichnan rapid change model which was investigated in 
Ref. [27|. The answer on this question can be treated as 
the first step on the way of investigating of the model 
with velocity field driven by the stochastic Navier-Stokes 
equation which is more complicated form mathematical 
point of view. The second aim is to analyze whether the 
finite correlation time of velocity field can lead to more 
complicated structure of critical dimensions than it was 
shown in Ref. (27j within the rapid-change model with 
small-scale anisotropy. 

The paper is organized as follows. In the first part of 
Sec.|lTl we give the precise formulation of used model. In 
the second part, we give the field theoretic formulation of 
the model and discuss corresponding diagrammatic tech- 
nique. In Sec. lIIIl we perform the ultraviolet (UV) renor- 
malization of the model, the renormalization constants 
are calculated in one-loop approximation, and the corre- 
sponding RG equations are derived. In Sec. lIVI we dis- 
cuss the stability of possible scaling regimes of the model 
which are governed by the corresponding infrared (IR) 
fixed points. In Sec.|Vl the renormalization of needed 
composite operators is done and their critical dimensions 
are found as functions of parameters of the model. Ob- 
tained results are reviewed and discussed in Sec.lVIl 



II. FIELD THEORETIC DESCRIPTION OF THE 
MODEL 

The advection of a passive scalar field 9(x) = 9(t, x) in 
an incompressible turbulent environment is described by 
the stochastic equation 

dtO + VidiO = u Q M + f, (1) 



where dt = d/dt, <9; = d/dxi, v$ is the molecular diffu- 
sivity coefficient (in what follows, a subscript denotes 
bare parameters of unrenormalized theory), A = d 2 is 
the Laplace operator, vi = Vi(x) is the z-th component of 
the divergence-free (owing to the incompressibility) ve- 
locity field v(x), and / = fix) is an artificial Gaussian 
random noise with zero mean and correlation function 

D e (x; x') = (f(x)f(x')) = S(t - t')C(r/L), 

(2) 

where the angular brackets (...) hereafter denote average 
over the corresponding statistical ensemble and L is an 
integral scale related to the stirring. The random noise is 
introduced to maintain the steady state of the system but 
the detailed form of the function C(r/L) in Eq. (J5J) will 
be inessential in our consideration. The only condition 
which must be satisfied by the function C(r/L) is that 
it must be finite and must decrease rapidly for r 3> L. 
In the problems related to the genuine turbulence the 
velocity field v(x) satisfies Navier-Stokes equation but, 
in what follows, we shall work with a simplified model 
where we suppose that the statistics of the velocity field 
is given in the form of a Gaussian distribution with zero 
mean and pair correlation function [3(J [H[ 

r fji 1 )rl d '\< 

Dl{x;x>) = (vi{x)vj(x')) = y^-^Ck) 

x ^(w,fc)e" l[w(t - t ' ) - k - (x - x ' )1 , (3) 

with 

D v (u,k) = -— . — - (4) 

(tOJ + UqVqK a r >){— VjJ + UqVqK a 1) 

where k = |k| and a transverse projector Py(k) reflects 
vectorial nature of the solenoidal velocity field. In the 
isotropic case it has the form of the simple transverse 
projector 

P y(k )=%-Ml. (5 ) 

In the anisotropic case the transverse projector becomes 
more complicated as it will be specified below (see also 
Ref. [13 )• In Eq- © D = g ^ is a positive ampli- 
tude factor and introduced parameter go plays the role 
of the coupling constant of the model. In addition, go 
is a formal small parameter of the ordinary perturba- 
tion theory. On the other hand, the parameter uq, in- 
troduced in the denominator of Eq. gives the ratio of 
turnover time of scalar field and velocity correlation time 
(see, e.g., Ref. 30] for details). The positive exponents 
e and r\ (e = 0(rj)) are small RG expansion parameters. 
Thus, we have a kind of double expansion model in the 
e — r\ plane around the origin e = r\ = 0. The coupling 
constant go and the exponent e control the behavior of 
the equal-time pair correlation function of velocity field 
(mean square velocity) or, equivalently, energy spectrum. 
On the other hand, the parameter uq and the second 
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exponent r\ are related to the frequency a; ~ u^vqU 2-11 
which characterizes the mode k (30l l49l. l50l. l5ll. |52| . Thus, 
in our notation, the value e — 4/3 corresponds to the 
celebrated Kolmogorov "two-thirds law" for the spatial 
statistics of the velocity field or, equivalently, "five-thirds 
law" for the energy spectrum, and rj — A/3 corresponds 
to the Kolmogorov frequency. Simple dimensional anal- 
ysis shows that go and uq, which we commonly term as 
charges, are related to the characteristic ultraviolet (UV) 
momentum scale A (or inner legth I ~ A -1 ) by the fol- 
lowing relations 
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A 



u Q 



A". 



(0) 



As was discussed in Introduction, in what follows, we 
shall take the velocity statistics to be anisotropic at all 
scales. For that purpose, we replace the ordinary trans- 
verse projector -FVj(k) in Eq. ^ with the general uni- 
axially anisotropic transverse tensor structure (see, e.g., 
Ref. [23 ): 

T tf (k) = a(ip)Pij(k) + &(VOPi.(k)n a ntP ti (k), (7) 

where n, is the i-th component of the unit vector n 
(n 2 = 1) which determines the distinguished direction 
of uniaxial anisotropy and ijj is the angle between the 
vectors k and n, so that n • k = k cos-0- It is well known 
that functions a(if>) and b(ip) can be decomposed into the 
d-dimensional generalization of the Legendre polynomi- 
als which are known as the Gegenbauer polynomials [53| , 
namely, 

00 00 
a (V0 = ^2 ai P 2 i {cos ip), b{ip) = ^2biP 2 i(cosip) 
1=0 1=0 

_ (8) 
(as was shown in Ref. [27[ the odd polynomials do not 
affect the scaling behavior). The necessary condition to 
have positively defined velocity correlator ([3]) leads to the 
following inequalities for these functions [271 ]: 



a(ip)>0, a(ip) + b(ip) sm 2 ip > 0. 



(9) 



But in practical calculations it is impossible to work with 
the general tensor structure as is defined in Eq. (J7J . The 
reason is, at least, because it contains infinite number of 
parameters a* and bj in the corresponding decomposition 
(|5J). Therefore, in what follows, we shall work with the 
simplest special case of the general uniaxial anisotropic 
transverse projector, namely, 

nk 



r i3 -(k) = [1 + ai 



Pij (k) + a 2 P is (k)n s n t P ti (k) , 

which is sufficient for investigation of principal properties 
of the uniaxial anisotropy (see the corresponding discus- 
sion in Ref. (27J ) . In this case, the inequalities © reduce 
into the requirements a\ > — l,a 2 > — 1- This special 
case represents nicely all main features of the general 
mod el (171 ). This can be seen from the analysis given in 
Ref. ^$7 



Let us briefly discuss two special limits of the con- 
sidered model ©, 0]) (see also Ref. They will be 
also studied in what follows. The first of them is the 
so-called rapid-change model limit when uo — > 00 and 
9o = 9o/ u o = const 



D v (u;,k)^g' v Q k- d - 2 ^, 



(11) 



and the second is the so-called quenched (time- 
independent or frozen) velocity field limit, which is de- 
fined by uq — > and g' ' = go/uo = const 



D v {w,k) -> g'^ 2 Trd(uj)k 



-d+2-2e 



(12) 



which is similar to the well-known models of random 
walks in a random environment with long-range corre- 
lations; see, e.g., Refs. IHlj. 

Using the well-known Martin-Siggia-Rose mechanism 
(see also, e.g., Refs. 0, HH) the stochastic problem 
(HI)-® can be treated as a field theory with action func- 
tional 



S(0,0',v) = -~ J dt 1 d d x 1 dt 2 d d x 2 



Vi (ti , Xi ) [£)?• (ti , Xi ; t 2 , x 2 )] Vj {t 2 , x 2 ) 
1 



1 j \ 



dti d a x x dt 2 <f x 2 
0'(ii,x 1 ) J D e (t 1 ,x 1 ;t 2 ,x 2 )0'(i 2 ,x 2 ) 
+ J dt d d x 9' [~d t - Vidi + v a A] 9 , (13) 

where 0' is an auxiliary scalar field, and D and D v are 
correlators (J2J) and ((3]), respectively. In action (fl3|) all 
required summations over the vector indices are under- 
stood. The second and the third integral in Eq. (fl3|) 
represent the DeDominicis-Janssen-type action for the 
stochastic problem (pTJ) , @ at fixed v, and the first inte- 
gral represents the Gaussian averaging over v. 

Model (fT3)) corresponds to a standard Feynman dia- 
grammatic technique with the bare propagators {99') 
and (viVj)o (in the time-momentum representation) 



{9{t,k)0'{t',k)) o - 
(w i (t,k)u j (t / ,k)) = 



8(t-t , )e-" ok2{t ~ t '\ (14) 
Do 



2u fc 1+2e 
x e-^^ofe 2 -" (t-f) Py .( k ) ; ( 15 ) 

where 9(t — t') is the step function, or (in the frequency- 
momentum representation) 



{9(uj,k)9'(u;,k)) 



1 



-iuj + v§k 2 ' 



(16) 



(viv^o = Tij(k)D v {uj,k), (17) 

where D v {uj, k) is given directly by Eq. ((4]). In the Feyn- 
man diagrams these propagators are represented by the 
lines which are shown in Fig.[T] (the end with a slash in 
the propagator {80') corresponds to the field 9', and the 
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(vv) Q = 
(0O'h = 




FIG. 1: (Left) Graphical representation of needed propaga- 
tors of the model. (Right) The triple (interaction) vertex of 
the model. Momentum k is flowing into the vertex via the 
auxiliary field 8' . 



end without a slash corresponds to the field 8) . The triple 
vertex (or interaction vertex) —O'vjdjO = 8'vjVj9, where 
Vj = ikj (in the momentum- frequency representation), is 
present in Fig.[TJ where momentum k is flowing into the 
vertex via the auxiliary field 8'. 

In the presence of anisotropy to have a multiplicatively 
renormalized model it is also necessary to introduce new 
counterterm of the form 8'(n-d) 2 8, which is absent in the 
unrenormalized action functional (| 13() . It means that the 
model given by action (| 1 3[) in its original formulation is 
not multiplicatively renormalizable, and in order to use 
the standard RG technique it is necessary to extend the 
model by adding the new contribution to the unrenor- 
malized action (1131). The extended action is 



S(9,8',v) = -\J dt l d d K 1 dt 2 d d *z 



Vi (ti , Xi ) [£)?• (t x , Xi ; t 2 , x 2 )] Vj (t 2 , x 2 ) 
1 



(18) 



dh cf xi dt 2 d a x 2 
0'(ii,x 1 )l? e (t 1 ,x 1 ;t 2 ,x 2 )0'(i 2 ,x 2 ) 
+ J dtd d x 9' [-d t ~ vA + v A + Xo ^o(n • d) 2 } 9 , 

Here xo is a new dimensionless unrenormalized parame- 
ter. The stability of the system implies the positivity of 
the total viscous contribution v$k 2 + xo^oink) 2 , which 
leads to the inequality Xo > ~~ 1- Its "physical" value 
is zero, but this fact does not hinder the use of the RG 
technique, in which it is first assumed to be arbitrary, 
and equality xo = is imposed as the initial condition in 
solving the equations for invariant variables. Below we 
shall see that the zero value of xo corresponds to certain 
nonzero value of its renormalized analog. 

For the action (fl8|) . the bare propagator in Eq. (fl6|) is 
replaced with 



?')o = 



1 



vok 2 + xo^o(nk) 



(19) 



The formulation of the problem through the action 
functional (fTgj) replaces the statistical averages of random 
quantities in the stochastic problem defined by Eqs. ([I]) 
and ([3]) with equivalent functional averages with weight 
expS*($), where $ = {8,8',v}. The generating func- 
tional of the total Green functions G(A) and connected 



TABLE I: Canonical dimensions of the fields and parameters 
of the model under consideration. 



F 


V 


e 


e' 


m, A, fi 




So 


u 


g,u,xo,x 


4 


-1 





d 


1 


-2 


2e + rj 


V 





d F 


1 


-1/2 


1/2 





1 











d f 


1 


-1 


d+1 


1 





2e + ri 


V 






Green functions W(A) are then defined by the functional 
integral 



G(A) = e w ^ = J e s (*)+^*, 



(20) 



where A(x) = {A e , A 6 , A v } represents a set of arbitrary 
sources for the set of fields $, = V8T)8'T)v denotes 
the measure of functional integration, and the linear form 
A<& is defined as 

A<P= { dx[A e {x)8{x)+A e ' (x)8'(x)+A v l (x)v l (x)}. (21) 



III. RENORMALIZATION GROUP ANALYSIS 

Using the standard analysis of canonical dimensions 
leads to the information about possible UV divergences 
in the model (see, e.g., Refs. [2(| HH). The dynamical 
model (Q~8]) belongs to the class of the so-called two-scale 
models [2lL [H, |24| . i.e., to the class of models for which 
the canonical dimension of some quantity F is given by 
two numbers, namely, the momentum dimension dp and 
the frequency dimension d F . To find the dimensions of all 
quantities it is convenient to use the standard normaliza- 



tion conditions d\ 

dt = 4 = 



-d* = i,d% = -d? = i,d% = d- 



0, and the requirement that each term of 
the action functional must be dimensionless separately 
with respect to the momentum and frequency dimen- 
sions. The total canonical dimension dp is then defined 
a,s dp — d F + 2d F (it is related to the fact that dt oc v^d 2 
in the free action (fT5|) with choice of zero canonical di- 
mension for i>o)- In the framework of the theory of renor- 
malization the total canonical dimension in dynamical 
models plays the same role as the momentum dimension 
does in static models. 

The canonical dimensions of our model are present in 
TablelH where also the canonical dimensions of the renor- 
malized parameters are shown. 

The necessity to work with the model based on the 
action (fT5)) instead of the action (JT3J) is given by the fol- 
lowing consideration. The model (|13p is logarithmic at 
e = r\ = (the coupling constants go and uq are dimen- 
sionless); therefore, in the framework of the minimal sub- 
straction (MS) scheme [2(| , which is always used in what 
follows, possible UV divergences in the correlation func- 
tions have the form of poles in e, 77, and their linear combi- 
nations. It is well known that the superficial divergences 
can be present only in the 1-irreducible Green functions 
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for which the corresponding total canonical dimensions 
are a nonnegative integer. Detail analysis of the possible 
divergences was done, e.g., in Ref. therefore we shall 
not repeat it here. This analysis shows that superficially 
divergent function of our model is only function (9'0)i—i r . 
From the action functional (|13|) one immediately obtains 
that the corresponding counterterms, which are needed 
to remove these divergences, must be proportional to two 
symbols d and, in the isotropic case, it is reduced to the 
structure 8' A9. However, in the anisotropic case, it is 
necessary to introduce possible anisotropic counterterm 
9'(n ■ d) 2 9 which is not present in the original action 
(p~3|) but which is generated during calculations. This is 
the reason why our starting action is the action given in 
Eq.flB). 

After this extension the model has become multiplica- 
tively renormalizable. It means that all divergences can 
be removed b y th e counterterms of the forms 9'A9 and 
9'(n-d) 2 9 [27I.T30I] . This can be explicitly expressed in the 
multiplicative renormalization of the parameters ga,uo, 
u , and Xo m the form 

v a = vZ v ,g a = g^ 2e+ ' n Z g , u a = uyPZ u , \o = X%x- (22) 

Here the dimensionless parameters g,u,v, and x are 
the renormalized counterparts of the corresponding bare 
ones, fj, is the renormalization mass (a scale setting pa- 
rameter), an artefact of the dimensional regularization. 
Quantities Zi — Zi{g, u, x; d; e, 77) are the so-called renor- 
malization constants and, in general, they contain poles 
in linear combinations of e and r\. 

The renormalized action functional has the following 
form: 

S R {6,0',v) = -\] dt 1 d d Xidt 2 d d x 2 

Vi (t 1 , xi ) [D^ (t 1 , xi ; t 2 , x 2 )] ~ 1 Vj (t 2 , x 2 ) 
+ \J dtid d x 1 dt 2 d d x 2 (23) 
, (< 1 ,x 1 )^ 9 (t 1 ,x 1 ;i 2 ,x 2 )0'(t 2 ,x 2 ) 
+ J dtd d x 9' [-d t - v l d l + vZ x A + xvZ 2 {n ■ d) 2 ] 9 , 

By comparison of the renormalized action (|23|) with defi- 
nitions of the renormalization constants Zi, i = g,u, v, x> 
which are given in Eqs. (|22p . we come to the relations 
among them: 

Z v = Z\, Z x = Z 2 Z 1 , Zg = Z 1 , Z u = Z-y , (24) 

The issue of interest is, in particular, the behavior of re- 
sponse functions, e.g., (6(x)6'(x')), correlation functions 
(9(xi)9(x 2 )...9(x n )}, and the equal-time structure func- 
tions 

SW(r) = ([0(i,x)-0(t,x')H, r=|x-x'| (25) 

in the inertial range specified by the inequalities I ~ 
1/A C r C L = 1/m (/ is an internal length). In the 



/ \ 

^e'e = — M h 



FIG. 2: The one-loop diagram that contribute to the self- 
energy operator £9/9. 

field theoretic formulation of our stochastic problem the 
angular brackets (...) mean functional average over fields 
8,0',v with weight exp(Sn). Independence of the origi- 
nal unrenormalized model of the scale-setting parameter 
fi of the renormalized model yields the RG differential 
equations for the renormalized correlation functions of 
the fields, e.g., 

PV+ ftdi-7^]<0(x,i)0(x',i')>H = O. (26) 

Here T> x = xd x stands for any variable x and the RG 
functions (the (3 and 7 functions) are given by the well 
known definitions [2(1 [5l[ ■ In our case, using the relations 
(|24p for the renormalization constants, they acquire the 
following form: 

7l = V„\nZi (27) 

for any renormalization constant Zi, and 

[3 g = V^g = g{-2e-J 1 + i 11 ), (28) 
(3 U = V ll u = u{-r } + ^i), (29) 
Px = v^x = xhi - 72). (30) 

The renormalization constants Z\ and Z 2 are deter- 
mined by the requirement that the one-particle irre- 
ducible Green function (0'9)i-i r must be UV finite when 
is written in the renormalized variables. In our case this 
means that it has no singularities in the limit e,i] — > 0. 
The one-particle irreducible Green function (9'6)\-i r is 
related to the self-energy operator T,g>g, which is ex- 
pressed via Feynman graphs, by the Dyson equation. In 
frequency-momentum representation it has the following 
form: 

(fl^i-ir = -iu + vop 2 + ^Xo(n-p) 2 -E e >e(u,p)- (31) 

Thus Z\ and Z 2 are found from the requirement that the 
UV divergences are canceled in Eq. (f3"Tj) after the substi- 
tution uq = vZ v , xo = x%x- This determines Z\ and 
Z 2 up to an UV finite contribution, which is fixed by the 
choice of the renormalization scheme. In the MS scheme 
all the renormalization constants have the form: 1 + 
poles in e, rj and their linear combinations. In one-loop 
approximation the self-energy operator Y,g:g is defined by 
Feynman diagram which is shown in Fig. [21 

It can be shown that in one-loop calculations it is 
enough to work with 77 = (see, e.g., Refs. [13, HH, [HJ for 
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details). This possibility essentially simplifies the evalu- 
ations of all quantities. Then the divergent part of the 
diagram given in Fig. [2] has only poles in e. Its explicit 
analytical form is given as follows (in rcnormalized pa- 
rameters and within one- loop approximation): 



Sd 



l l 



(27r) d 2u{\ + u) d(d + 2)e 
x [p 2 A+(n-p) 2 B] , 



with 



(32) 



The fixed point of the RG equations is defined by (3- 
functions, namely, by requirement of their vanishing. In 
our model the coordinates g* , it* , \* of all possible fixed 
points arc found from the system of three equations 



(3 g {g*,u*,x*) = P u (g*,u*,x*) = /3 x (g*,u*,x*) = 0. 



(39) 



The /3-f unctions g , U , and j3 x are defined in Eqs. 
(|29)l . and ((30|) . To investigate the IR stability of a fixed 
point it is enough to analyze the eigenvalues of the matrix 
of the first derivatives: 



(33) 



A = (l + a 1 )d(d + 2) 2 F 1 (l;li;-^- 

\ 2 2 1 + U 

+ (a 2 -a 1 d-l)(d+2) 2 F 1 (l;\;l + ^-^- 

\ 2 2 1 + u 

+ ( ai -a 2 )(d+l) 2 F 1 (v,l-,2+~-,^- 
\ 2 2 1 + u 

B = _(l + ai )d(d + 2) 2 Fifl;^;-=^- 

\ 2 2 1 + u 

-{axil - 2d) + a 2 -d](d+ 2) 2 Fx ( 1; \;l + % — — 

\ 2 2 1 + It 

-(«! - a 2 )d(d + 1) 2 F, (l; 1; 2 + |; . (34) 

where = 2ir d / 2 /T(d/2) denotes surface of the <i- 
dimensional unit sphere and 2 Fi(a,b, c, z) = 1 + + 

o(a+l)6(6+l) „2 



c(c+l)-l-2 



z + . . . represents the corresponding hyper- 
geometric function. 

In the end, the renormalization constants Z\ and Z 2 
are given as follows 



Z x = 1 
Z 2 = 1 



2u(l + u) d(d + 2) e' 

g 1 g 

2u(l + u) d(d + 2)x e ' 



(35) 
(36) 



where we have introduced new notation g = gSd/ {2ii) d . 

Now using the definition of the anomalous dimensions 
71 and 72 in Eq. ([27f one comes to the following expres- 
sions: 



7i 
72 



9 



2u(l + u)d(d + 2) 
9 



2u{l + u)d{d + 2)x 



A. 
B. 



(37) 
(38) 



In the next section we shall use these results for investi- 
gation of possible scaling regimes of the model. 



IV. FIXED POINTS AND SCALING REGIMES 

Possible scaling regimes of a renormalized model are 
directly given by the infrared (IR) stable fixed points of 
the corresponding system of the RG equations [2(1 Hl[ • 



dP g /dg df3 g /du dp g /d X 
fiij = I df3 u /dg df3 u /du d/3 u /dx 
8f3 x /dg df3 x /du d(3 x /d X 



(40) 



The IR asymptotic behavior is governed by the IR stable 
fixed points, i.e., those for which real parts of all eigen- 
values are nonnegative. 

First of all, we shall study the rapid-change model 
limit: u — > 00. In this regime, it is convenient to make 
transformation to new variables, namely, w = 1/it, and 
g' = g/u 2 [30I ] . with the corresponding changes in the (3 
functions: 



0g. = g'{-2e + 77 + 7l ), 
Pw = w(t]-^x), 



(41) 
(42) 



while (3 X is unchanged, i.e., it is given by Eq. (|30p . In this 
notation the anomalous dimensions 71 and 72 acquire the 
following form: 



7i 

72 



2(l + w)d(d + 2) 
-/ 



2(1 + w)d(d + 2) X ' 



(43) 
(44) 



where again g' = g'Sd/(2ir) d , and A' and B' acquire the 
form 

A' = (1 + a x )d{d + 2) 2 Fx (l; \; ~; ^) (45) 

+ (a 2 ~ axd - l)(d + 2) 2 Fx fl; ~; 1 + ~; ^ 

\ 2 2 1 + w 

+ (ax-a 2 )(d+l) 2 Fx (l-,\-2+~-^- 
\ 2 2 1 + w 

B> = -(1 + ax)d[d + 2) 2 Fx ^) (46) 

-[ai(l - 2d) + a 2 -d](d + 2) 2 F X fl; 1 + J -=£H1 

\ 2 2 1 + w 

-(ai-a a )(d + l)<i 3 fi (l; I; 2 + ^; 

V 2 2 1 + tu 



In the rapid-change model limit w — > (it — > 00) we 
are coming to the result of Refs. [27[ with the anomalous 
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dimensions 71 and 72 of the form 



71 = lim 



w->o 2(1 + w)d(d + 2) 
9' 



A' 



(47) 



72 



2d(d + 2) 
lim 



[(d-l)(d + 2)+a 1 (d+l) + a 3 ] 1 



««-o 2(1 + w)d(d + 2)x 



D' 



2d(d + 2) x 



[-2 ai + (d 2 -2)a 2 ]. 



(48) 



For completeness we shall briefly discuss this spacial case. 
In this limit we have two fixed points denoted as FPI and 
FPU. The first fixed point is trivial, namely 



FPI 



= 9* = 0, 



(49) 



with arbitrary \* and 7i = 0, 72 =0. The corresponding 
"stability matrix" is triangular with diagonal elements 
(eigenvalues) : 



Ai = -2e + 77, A 2 = t], A 3 = 0. 



(50) 



The region of the IR stability is shown in Fig. [31 The 
second point is defined as 



FPU 



0, 



2d(d + 2)(2e-r/) 



X* 



(d+2)(d-l)+ai(d + l) + a 2 

-2a! + a 2 (d 2 - 2) 
(d + 2)(d-l) + ai(d+l) + a 2 



(51) 
,(52) 

.(53) 



with 7^ = 72 = 2e — 77. The triangular matrix £1 has the 
following eigenvalues (diagonal elements) 

Ai = 2e - n, A 2 = 2e - 77, A 3 = -2e + 2<q. (54) 

The region of the IR stability of this fixed point is shown 
in Fig.d 

Now let us analyze the "frozen regime" with frozen 
velocity field. It is mathematically obtained from the 
model under consideration in the limit u — > 0. To study 
this transition it is appropriate to change the variable g 
to the new variable g" = g/u [3(J. Then the /3 9 function 
is transformed to the following one: 



/V= ff /, (-2e + 2 7 i), 



(55) 



while (3 U and (3 X functions are not changed, i.e., they are 
the same as the initial ones given by Eqs. |29)) and (f30j) . 
In this notation the anomalous dimensions 71 and 72 have 
the form 



71 
72 



2(1 + u)d(d + 2) 
9" 

2(l + u)d(d + 2) X 



A. 



(56) 
(57) 



5 

T|4 

3 
2 
1 

-1 
-2 
-3 
-4 
-5 



FPI 


£=r|/2/ FPII/ 
/ / e=r| 

/ /^ FPV 


FPIII 


FPIV 



5 -4 -3 -2 -1 1 2 3 4 5 



FIG. 3: The "phase" diagram of the fixed points of the model 
(see the text for details). 



where, as always, g" — g" Sd/ (2ir) d . In the limit u — » 
the anomalous dimensions 71 and 72 acquire the following 
form: 



71 
72 
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2d(d + 2 
9" 



■A" 



B" . 



2d(d + 2) X 
where A" and B" are given as follows 

A" = (l + a 1 )d(d + 2) 2 F 1 (l-ht- x 



(58) 
(59) 



B" 



+ (a 2 - aid - l)(d + 2) 2j Fi ^1; -; 1 + - ; ~x 

+(ai-a 2 )(d+l) 2 Fi fl-L^+^-x") , (60) 

1 d 



-(l + Ql )d(d+2) 2 Fa ^l;-;-;- X 
-[ai(l - 2d) + U2 - d]{d + 2) 2 fi ^1; i; 1 + ^; - X 
-(ai-a 2 )(d+l)d 2 Ji A;i;2+^;- x V(61) 



The system of /3 functions (|29|). p0|) . and (J55j) exhibits 
two fixed points, denoted as FPIII and FPIV, related to 
the corresponding two scaling regimes. One of them is 
again trivial, namely, 



FPIII 



= g'l = 0, 



(62) 
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with arbitrary x* and 7* = 72 — 0- The eigenvalues of 
the corresponding matrix O are 

Ai = -2e, A 2 = -t?, A 3 = 0. (63) 

Thus this regime is IR stable only if both parameters s, 
and rj are negative simultaneously as can be seen in FigOJ 
The second, non-trivial, point is 

FPIV: = 0, (64) 

where A'l is A" given in Eq. (|60p taken at fixed point, 
i.e., x is replaced by x* which is given only implicitly by 
the equation 

X*K - K = 0, (66) 

where B" is B" given in Eq. (|61[) taken at the fixed point. 

Straightforward analysis shows that to have g" > 
together with x» > — 1 one must suppose e > 0. It is 
the only condition related to the coordinates of the fixed 
point. 

The IR stability of the fixed point is again given by 
the f2 matrix, namely, by the positive values of real parts 
of its eigenvalues. It is triangular in this case, thus its 
eigenvalues are given directly by the diagonal elements. 
The eigenvalues are 

A a - 2e, (67) 
A 2 = e-V, (68) 

(69 » 

Here A3 has rather complicated explicit form but it can 
be numerically shown that A3 is always positive for 0:1,2 > 
— 1, e > 0, and d > 0. The region of stability of this fixed 
point is shown in Fig.[3J 

Now let us turn to the most interesting scaling regime 
with finite value of the fixed point for the variable u. By 
short analysis one immediately concludes that the system 
of equations 



g = g(-2e - 77 + 3 7l ) = 0, (70) 
Pu = u(-77 + 7i) = 0, (71) 
Px = X(n - 72) = 0. (72) 

can be fulfilled simultaneously for finite values of g and 
u only when the parameter e is equal to 77: e — rj. In 



this case the function (3 g is proportional to function {3 U . 
As a result we have not one fixed point but a set of fixed 
points g*, x* that depend on arbitrary parameter > 0. 
The value of the fixed point for the variable g in one-loop 
approximation is given as follows (we denote it as FPV) 

FPV : * = Ml + *.)d{d + 2)e ^ 



where A* is A from Eq. (|33f with u and x replaced 
by it, and x*i respectively On the other hand, x* is 
again known only implicitly and it can be obtained from 
Eq. (|T2"|) which is equivalent to the condition 

7i= 7a, (74) 

where 7^ , 72 are 71 , 72 given by Eqs. (|3"7|) and (|38[) where 
g, u are replaced by g* and u* respectively. 

The eigenvalues of the corresponding stability matrix 
are 

Ai = (75) 
A 2 , 3 = \[C± VC 2 -4D], (76) 

where 

C = 3e + x*<9 x (7i ~ 72)1* + 7i*<9„7i|» , 
D = 3ex*<9 x (7i -72)!* 

-X*u*(du7i <9 X 72 - <9 x 7i <9„7 2 )|* 

where |* means that the quantity must be taken at the 
fixed point. It can be shown numerically that for any 
positive values of and for all possible values of the 
anisotropy parameters 01,2 the eigenvalues A2 and A3 are 
always greater then zero. Therefore, the corresponding 
fixed point is IR stable and satisfy stability condition. It 
corresponds to the line e = rj in Fig.[3l where the regions 
of stability for all possible fixed points are shown. 

As was already mentioned (see the previous section) 
the issue of interest are especially multiplicatively renor- 
malizable equal-time two-point quantities G(r) (see also, 
e.g., Ref. [3fJ|). Examples of such quantities are the equal- 
time structure functions in the incrtial interval as they 
were defined in Eq. I]25p. The IR scaling behavior of the 
function G(r) (for r//> 1 and any fixed r/L) 

G{r) ~v^ G r dG (r/l)- AG R{r/L) (77) 

is related to the existence of IR stable fixed points of the 
RG equations (see above) . In Eq. (|77|) and do are cor- 
responding canonical dimensions of the function G (the 
canonical dimensions of the model are given in Scc. lIII[) . 
R(r/L) is the so-called scaling function, which cannot be 
determined by the RG equation (see, e.g., Ref. [2l|]). and 
Ag is the critical dimension defined as 

A G = 4 + A w dg + 7£. (78) 

Here is the fixed point value of the anomalous dimen- 
sion 7g = fid^hiZQ, where Zq is the renormalization 
constant of the multiplicatively renormalizablc quantity 
G, i.e., G = Z G G R and A w = 2 - 7* is the crit- 
ical dimension of the frequency with 7* = 7^ which is 
defined in Eq. (|3"T|) and ^{ means that is taken at the 
corresponding fixed point. From above discussion of the 
possible scaling regimes we have 

{ 2e-r] for FPU \ 
7* ee i = \ e for FPIV \ . (79) 
[ e = rj for FPV J 
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r(i) 




FIG. 4: Graphical representation of the one-loop correction 
to Tjv in Eq. (1M[) . 



We are working only in one-loop approximation but the 
anomalous dimension 7 * is already exact for all fixed 
points at one-loop level [30l. |34| . i.e., it has no loop cor- 
rections of higher order, therefore the critical dimensions 
of frequency u) and of fields <I> = {v, 9, 9' } ar e also found 
exactly at one- loop level approximation [30]. In our no- 
tation they read 



and 



1-7:, A e 



2 - 2e + 77 for FPU 
2 - e for FPIV 
2 - e = 2 - 77 for FPV 



(80) 



-l+ 7 :/2, A e , =d+l- 7 ;/2. (81) 



The renormalized function G fl must satisfy the RG 
equation of the form 

(V RG + lG )G R (r) = 0, (82) 

with operator T>rg given explicitly in Eq. (|26[) , namely, 



V 



RG 



(83) 



The difference between the functions G and G R is only in 
the normalization, choice of parameters (bare or renor- 
malized), and related to this choice the form of the per- 
turbation theory (in go or in <?). The existence of a non- 
trivial IR stable fixed point means that in the IR asymp- 
totic region r/l > 1 and any fixed r/L the function G(r) 
takes on the self-similar form given in Eq. ([77)) . As was 
already mentioned the scaling function R(r/L) is not de- 
termined by the RG equation itself. The dependence of 
the scaling functions on the argument r/L in the region 
r/L <C 1 can be studied using the well-known Wilson 
operator product expansion (OPE) [H HH, [H, It 
shows that, in the limit r/L — > 0, the function R{r/L) 
can be written in the following asymptotic form: 



R{r/L) = Y,C Fi (r/L) (r/L) 1 



(84) 



where C Fi are coefficients regular in r/L. In general, 
the summation is implied over certain renormalized com- 
posite operators Fi with critical dimensions A. Fi . In 



A[2,2R a =0, d=3 




a. 



FIG. 5: Dependence of the critical dimension A[2, 2]/f on 
anisotropy parameter ot\ (02 = 0) for different fixed point 
values of the parameter u: u* — (frozen limit) - solid line, 
u* — 0.5 - dash line, u* — 1 - dot line, u* — 5 - dash dot line, 
u* = 00 (rapid-change model limit) - dash dot dot line. The 
small figure shows details that are not visible in basic figure. 



the case under consideration the leading contribution 
is given by operators Fi having the form F[N,p] — 
9^0 ■ ■ ■ d lp 9(di9di9) n with N =p + 2n. In the next sec- 
tion we shall consider them in detail, where the complete 
one-loop calculation of the critical dimensions of the com- 
posite operators Fn will be presented for arbitrary values 
of N, d 7 u, and 0:12- 



V. CRITICAL DIMENSIONS OF COMPOSITE 
OPERATORS AND ANOMALOUS SCALING 

A. Operator product expansion 

According to the OPE [H H [H Si], the equal- 
time product Fi (x') F2(x") of two renormalized com- 
posite operators [53] at x = (x' + x")/2 = const and 
r = x' — x" — > can be written in the following form: 



Fx{x')F 2 (x" 



C Fi {v)Fi{^t), 



(85) 



where the summation is taken over all possible renor- 
malized local composite operators Fi allowed by sym- 
metry with definite critical dimensions A p t , and the 
functions C Fi are the corresponding Wilson coefficients 
regular in L~ 2 . The renormalized correlation function 
(Fi(x')F 2 (x")) can now be found by averaging Eq. ([55]) 
with the weight expS'^ with S R from Eq. The 
quantities (Fi) appear on the right-hand side and their 
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A[2,2]/t a =0, d=3 




A[3,1R a=0, d=3 




a 



FIG. 6: Dependence of the critical dimension A[2, 2]/f on 
anisotropy parameter 012 («i = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



FIG. 8: Dependence of the critical dimension A[3, on 
anisotropy parameter ai (02 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



A[2,2]/4, a 1= a 2 , d=3 



A[3,1R a 1= 0, d=3 





FIG. 7: Dependence of the critical dimension A[2, 2]/f on 
anisotropy parameter at\ = 012 for different fixed point values 
of the parameter u (for notation see the caption in Fig.©. 



FIG. 9: Dependence of the critical dimension A[3, on 
anisotropy parameter 012 (ai = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



asymptotic behavior in the limit L 1 — > is then found 
from the corresponding RG equations and has the form 
(Fi) oc L~ Ap i . 

From the OPE one can find that the scaling func- 
tion R(r/L) in the representation (|77p for the correla- 
tion function Fi(x')F2(x") has the form given in Eq. 



where the coefficients Cf, are regular in (r/L) 2 . 

It is well known that the specific feature of the turbu- 
lence models is the existence of operators with negative 
critical dimensions (the so-called " dangerous" operators) 
[U, HI HE HI- Their presence in the OPE deter- 
mines the IR behavior of the scaling functions and leads 



A[3,1R a=a 2 , 6=3 

0.5-j — ' ' 1 ' ' ■ 1 

0.14H ■ ' ■ ' ■ ' — 




FIG. 10: Dependence of the critical dimension A[3, l]/f on 
anisotropy parameter ot\ = a? for different fixed point values 
of the parameter u (for notation see the caption in Fig. [J}. 



A[3,3R oc 2 =0, d=3 
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FIG. 11: Dependence of the critical dimension A[3, 3]/£ on 
anisotropy parameter ai (a.2 — 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 
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A[3,3]/S, a 1= 0, d=3 

26 i , 1 , 1 , 1 




1.7-L-, , , , , , , , V 

5 10 15 20 



FIG. 12: Dependence of the critical dimension A[3, 3]/£ on 
anisotropy parameter 012 (on = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



A[3,3]/^, a =a , d=3 




5 _, 10 15 20 



FIG. 13: Dependence of the critical dimension A[3, 3]/£ on 
anisotropy parameter ot\ — a-z for different fixed point values 
of the parameter u (for notation see the caption in Fig. [5}. 



to their singular dependence on L when r/L — > 0. At 
this point the turbulence models are crucially different 
from the models of critical phenomena, where the lead- 
ing contribution to the representation (|77p is given by the 
simplest operator F = 1 with the dimension A^? = 0, and 
the other operators determine only the corrections that 



vanish for r/L — > 0. If the spectrum of the dimensions 
Ap i for a given scaling function is bounded from below, 
the leading term of its behavior for r/L — > is give n by 
the minimal dimension. As was discussed in Ref. [3fJ, the 
model under consideration belongs to this case for small 
enough values of the exponents e, 77. 
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A[4,0R a =0, d=3 
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FIG. 14: Dependence of the critical dimension A[4,0]/£ on 
anisotropy parameter ai (ct2 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



FIG. 15: Dependence of the critical dimension A[4,0]/£ on 
anisotropy parameter 012 (02 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



In what follows, we shall concentrate on the equal- 
time structure functions of the scalar field as defined in 
Eq. (I25p . The representation ([77)) is valid with the di- 
mensions = -N/2, d G = -N, and A G = NA e = 
N(—l + 7^/2). In general, not only do the operators 
which are present in the corresponding Taylor expansion 
are entering into the OPE but also all possible operators 
that admix to them in renormalization. In the present 
anisotropic model the leading contribution of the Tay- 
lor expansion for the structure functions (|25[) is given by 
the tensor composite operators constructed solely of the 
scalar gradients 



F[N lP ] = d ll e---d lp 6{d % ed l o) n , 



(86) 



where N = p + 2n is the total number of the fields 9 
entering into the operator and p is the number of the 
free vector indices (see, e.g., Ref. [27lj ) . 



A[4,0R a 1= a 2 , d=3 




B. Composite operators F[N,p]: renormalization 
and critical dimensions 

Let us briefly discuss renormalization of the composite 
operators (|56"|) . A complete and detailed discussion of the 
renormalization of the composite operators is given in 
Ref. [26| . Therefore, we shall discuss only basic moments 
necessary to present explicit expressions for composite 
operators. 

The necessity of additional renormalization of the com- 
posite operators (|86"|) is related to the fact that the 
coincidence of the field arguments in Green functions 



FIG. 16: Dependence of the critical dimension A[4,0]/£ on 
anisotropy parameter ot\ — Q2 for different fixed point values 
of the parameter u (for notation see the caption in Fig. [5}. 



containing them leads to additional UV divergences. 
These divergences must be removed by special kind of 
renormalizati on p rocedure which can be found, e.g., in 
Refs. [HI, l2l| . where their renormalization is stud- 
ied in general. As for the renormalization of composite 
operators in the models of turbulence it is discussed in 
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A[4,2]/£, cc^O, d=3 
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FIG. 17: Dependence of the critical dimension A[4, 2]/£ on 
anisotropy parameter ai (012 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



A[4,2R a 2 =0, d=3 




FIG. 18: Dependence of the critical dimension A[4, 2]/£ on 
anisotropy parameter 012 (a.2 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 

Refs. [H, EH- Besides, typically, the composite opera- 
tors are mixed under renormalization. Therefore, let us 
briefly discuss this issue [2TI ] . 

Let F = {F a } be a closed set of composite opera- 
tors which are mixed only with each other in renormal- 
ization. Then the renormalization matrix Zp = {Z a p} 



A[4,2]/t a=a 2 , 6=3 

i i i i i i i 

1.0- I 



0.412- \ 
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FIG. 19: Dependence of the critical dimension A[4, 2]/£ on 
anisotropy parameter ot\ = 02 for different fixed point values 
of the parameter u (for notation see the caption in Fig.[SJ. 

and the matrix of corresponding anomalous dimensions 
7j? = {jap} for this set are given as follows 

K = Y1 Z *P F *' 7f = Z^D^Zp. (87) 

P 

Renormalized composite operators are subject to the fol- 
lowing RG differential equations 

+ m-luV u )F^ = -J2l^F^ (88) 

which lead to the following matrix of critical dimensions 
A F = {A al3 } 

A F = d F + A u d F +~/ F , A u = 2- 7 ;, (89) 

where d k F a d F are diagonal matrices of corresponding- 
canonical dimensions and "f F is the matrix of anomalous 
dimensions (|57|) taken at the fixed point. In the end, the 
critical dimensions of the set of operators F = {F a } are 
given by the eigenvalues of the matrix Ap. The so-called 
" basis" operators that possess definite critical dimensions 
have the form 

Fl as =£EV^ , (90) 


where the matrix Up = {U a p} is such that A' F = 
UpApUp 1 is diagonal. 

As was already mentioned, in what follows, the cen- 
tral role is played by the tensor composite operators 
di t 6 ■ ■ ■ di 6 (pi0di8) n , constructed solely of the scalar 
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A[5,1]/£, a 2 =0, d=3 A[5,1]/£, a=a 2 , d=3 




FIG. 20: Dependence of the critical dimension A [5, on 
anisotropy parameter ai (ct2 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



A[5,1]/4, a 1= 0, d=3 

-0.6 T—i ■ ' ■ 1 ■ ' 




5 10 15 20 

FIG. 21: Dependence of the critical dimension A[5, l]/f on 
anisotropy parameter «2 (o?i = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 

gradients. It is convenient to deal with the scalar op- 
erators obtained by contracting the tensors with the ap- 
propriate number of the uniaxial anisotropy vectors n 

F[N,p] = [(n-d)6] p (di6di6) n , N = 2n + p. (91) 



FIG. 22: Dependence of the critical dimension A[5, l]/f on 
anisotropy parameter ot\ = oli for different fixed point values 
of the parameter u (for notation see the caption in Fig.[SJ. 

Detail analysis shows that the composite operators (j9"Tj) 
with different N are not mixed in renormalization, 
and therefore the corresponding renormalization matrix 
Z[N, P ][N'y] is m fact block-diagonal, i.e., Z[ N , P ][N', P '] = 
for N' ^ ' N [13. 

In the isotropic case, as well as in the case when large- 
scale anisotropy is present, the elements Z\n,p] \N,p'] van- 
ish for p < p', thus the block Zmy my] is in fact triangu- 
lar along with the corresponding blocks of the matrices 
Uf and A F from Eqs. ([90]) and ([89]) . In the isotropic 
case it can be diagonalized by changing to irreducible 
operators (scalars, vectors, and traceless tensors), but 
even for nonzero imposed gradient its eigenvalues are the 
same as in the isotropic case. Therefore, the inclusion of 
large-scale anisotropy does not affect critical dimensions 
of the operators (ji?Tj) . On the other hand, in the case of 
small-scale anisotropy, the operators with different val- 
ues of p mix heavily in renormalization, and the matrix 
Z[N,p] [N,p'] is neither diagonal nor triangular here and 
one can write 

LAT/2J 

F[N,p] = Y, Z \n,p] [N,N-2i]F R [N, N-21], (92) 

1=0 

where [N/2\ means the integer part of the N/2. There- 
fore, each block of renormalization constants with given 
N is an ([N/2\ + 1) x {[N/2\ + 1) matrix. Of course, 
the matrix of critical dimensions (|89[) . whose eigenval- 
ues at IR stable fixed point are the critical dimensions 
A[N,p] of the set of operators F[N,p], has also dimen- 
sion (L-/V/2J + 1) x ([N/2\ + 1). 

Now let us turn to the calculation of the renormaliza- 
tion constants ^[jv,p] [JV,p'l in the one- loop approximation 
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A[5,3R a =0, d=3 A[5,3R a=a 2 , 6=3 
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FIG. 23: Dependence of the critical dimension A [5, 3]/£ on 
anisotropy parameter ai (012 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



A[5,3]/^, a 1= 0, d=3 




FIG. 24: Dependence of the critical dimension A[5, 3]/£ on 
anisotropy parameter «2 (o?i = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



in our model. We shall proceed as in Refs. [23, [30] • Let 
T(x;9) be the generating functional of the 1-irreducible 
Green functions with one composite operator F[N,p] 
from Eq. ([9"T]) and any number of fields 9. We shall be 
interested in the iV-th term of the expansion of T(x; 6) 



FIG. 25: Dependence of the critical dimension A[5, 3]/£ on 
anisotropy parameter ot\ = 02 for different fixed point values 
of the parameter u (for notation see the caption in Fig.[SJ. 



in 9, which we denote T^(x; 9); it has the form 

T N (x;9) = i J dxi--- j dx N 9{xi)---9(x N ) 

x(F[N,p](x)9(x 1 )---9(x N )) 1 _ il , (93) 
and in the one-loop approximation it is given as 

T N = F[N,p]+T (1) , (94) 

where I^ 1 ) is given by the analytical calculation of the 
diagram in Fig. 01 and the first term in Eq. ([9~4"|) represents 
"tree" approximation (see also Ref. [271]). 

The black circle with two attached lines in the diagram 
in Fig.[4]denotes the variational derivative V(x; xi,X2) = 
S 2 F[N,p]/S9(xi)S9(x2)- It can be represented in the fol- 
lowing convenient form (2?| 



V{x; xi,x 2 ) 



diS(x — x\) djS(x — X2) 
d 2 



dajdcij 



(na) p (a 2 ) n 



(95) 



where a constant vector <Zj will be substituted with &i9(x) 
after the differentiation. Analytical form of the the dia- 
gram in Fig. [4] (without the symmetry factor 1/2) is the 
following: 

dxi ■ ■ ■ J dx A V(x; x\, x 2 )(9(xi)9 / (x 3 )} 

x(9(x 2 )9'(x i )) {v k (x 3 )v l (x i )) d k 9(x 3 )d l 9(x i ), (96) 

where the bare propagators are given in Eqs. ([3]), (fl3|) and 
the derivatives are related to the ordinary vertex factors 
shown in Fig.[TJ 
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A[6,0]/£, a 2 =0, d=3 A[6,0]/^ a=a 2 , 6=3 
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FIG. 26: Dependence of the critical dimension A [6, 0]/£ on 
anisotropy parameter ai (012 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



FIG. 28: Dependence of the critical dimension A[6,0]/£ on 
anisotropy parameter ot\ = 02 for different fixed point values 
of the parameter u (for notation see the caption in Fig.[SJ. 



A[6,0]/5, a 1= 0, d=3 




FIG. 27: Dependence of the critical dimension A[6,0]/£ on 
anisotropy parameter 012 (ai = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



We are interested in the UV divergent part of the 
expression (|96p which is needed for determination of 
the corresponding renormalization constants. But the 
needed UV divergent part is proportional to the poly- 
nomial built of N gradients di6{x) at a single spacetime 
point x, and all of them have been already extracted from 



the |96|) , namely, N — 2 gradients are given by the vertex 
(I95p and the others 2 gradients are given by the ordinary 
vertex factors in Fig.[TJ This important point from the 
view of calculations allows us to replace the gradients 
with the constant vectors a. Therefore, in the end, the 
divergent part of expression (|96p can be written in the 
following compact form: 



d 2 



audi 



da.tdaj 



(na)> 2 )' 



X 



ij, kl 1 



(97) 



with 



X 



ij, kl 



dx 3 J dxi d l (6(x)6'(x 3 )}a 
xd j (9(x)9'(x i )} (v k (x 3 )vi( X4l )) , (98) 



or, in the momentum-frequency representation (suitable 
for the further calculations), after integration over the 
frequency, 



X 



ij, kl 



dk 



[2ttY (k 2 + m 2 ) d / 2 +> 
1 1 



-T kl (k) 



A: 2 +x(nk) 2 fc 2 (l + u) + x(nk) 



,(99) 



with Dq from Eq. ^ and Tm from Eq. (|10p (we again 
use the possibility to work with i] = within one-loop 
approximation [30, [3l|). Expression ([99]) can be decom- 
posed into some tensor structures (see, e.g., Ref. and 
after rather long but direct calculations we are coming to 
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A[6,2R a 2 =0, d=3 
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FIG. 29: Dependence of the critical dimension A [6, 2]/£ on 
anisotropy parameter ai (ct2 = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



A[6,2]/^, cc^O, d=3 
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FIG. 30: Dependence of the critical dimension A[6, 2]/£ on 
anisotropy parameter 012 (ai = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



the following result for the quantity denned in Eq. (|97l) : 

A (-T -iQi F[N,P- 2] + O2 F[N,p] 
{2Tr) a 4ir Vm/ e 

+Q 3 F[N,p + 2} + Q A F[N,p + 4}}, (100) 

where we have substituted the unrenormalized quantities 
with the renormalized one, <Zj have been replaced with 
the gradients di9{x) (thus they again form the operators 
F[N, q], with q = p — 2,p,p + 2,p + 4), and the following 
notation was applied for the corresponding coefficients 



3=0 



H, - 



1 + u 



Gj), »=1,...,4, (101) 



where Hj and Gj are the hypergeometric functions of the 
following form: 



H 4 



,1 d 



2 

d 



2F1 [ 2 



x 



l + u 



with j = 0, ...,3, and coefficients for i = 1, ...,4 and 
j = 0, 3 are given in Appendix A. 

Using the standard renormalization procedure the 
renormalization constants Z\jf p] [N,p>] defined in Eq. (|92[) 
are found from the requirement that function (|94p is UV 
finite (contains no poles in e) when is written in renor- 
malized variables and with the replacement F[N,p] — > 



A[6,2]/^ a=a 2 , 6=3 




FIG. 31: Dependence of the critical dimension A[6, 2]/£ on 
anisotropy parameter ot\ — Q2 for different fixed point values 
of the parameter u (for notation see the caption in Fig.[5j. 
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F R [N,p}. In the end, from Eqs. (J94J) and (flOPf we have 



J [N,p][N,p-2] 



Z[N,p][N,p] — 1 



z 



[N,p][N,p+2] 



'[N,p][N,p+4] 



fJ 



8u 2 s 

9 
8u 2 e 



(102) 
(103) 
(104) 
(105) 



with coefficients Qi given in Eq. (jlOip . Using the defini- 
tion of the matrix of anomalous dimensions 7[jv,p] [JV'.p'l 
given in Eq. (|87[) we are coming to the following result 



Qu 





g_ 


7[N,p][N,p-2] — 


Au 2 




9_ 


l[N,p][N,p] = 


Au 2 




9_ 


7[N,p][N,p+2] = 


Au 2 


l[N,p][N,p+4] = 


9_ 


Au 2 



Qu 



and the desired matrix of critical dimensions 
form 



(106) 



has the 



A[JV, P ][iV )P '] = ^7^/2 + 7*iv,p][iv,p']= ( 107 ) 

where the asterisk means that the quantities are taken 
at the corresponding fixed point (see Sec.lIV|) and 7* is 
given in Eq. (|79[) . The nonzero one- loop contribution to 
the matrix of critical dimension (|107[) is represented by 
Eqs.fUnil) with Qi,i = 1,...,4 defined in Eq. CEOI])- It 
means that the matrix elements of the matrix 7[at jP ] [N'y] 
other than given in Eq. (|106p are equal to zero. It can be 
seen immediately that the matrix of critical dimensions 
depends on the anisotropy parameters a\ and 012 and, 
what is now more interesting and important here, on the 
parameter u (see below). 

In the end, the critical dimensions A[N,p] are given 
by the eigenvalues of the matrix (|107[) . The simplest sit- 
uation occurs in the isotropic limit with a\ = 012 = 
and, correspondingly, \* — 0- In this case, one comes to 
the triangular matrix, therefore its eigenvalues are given 
directly by the diagonal elements. But more interesting 
is the fact that within the isotropic model we have the 
same eigenvalues of the matrix of critical dimensions for 
all fixed point values of u* , i.e., the eigenvalues are inde- 
pendent of u at the fixed point, namely, 



p(p-l) -n(d-l)(d + N + p) 



(d-l)(d + 2) 

r* (108) 

where £ is given in Eq. (|79|) (see, e.g., Ref. [30j for details). 
As a result, it means that within the one-loop approxi- 
mation there is no difference between general model with 
finite time correlations and its two special limits, namely, 
Kraichnan's rapid change limit and the frozen limit of the 
model as for the anomalous behavior of the equal-time 



A[7,m, a=0,d=3 




FIG. 32: Dependence of the critical dimension A[7, l]/£ on 
anisotropy parameter ai (02 — 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



A[7,1]/4, a 1= 0, d=3 




FIG. 33: Dependence of the critical dimension A[7, on 
anisotropy parameter 012 (ai — 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



structure functions (it, of course, also holds for the other 
equal-time correlation functions). 

The situation is different when presence of small scale 
anisotropy is supposed. In this case, the matrix of critical 
dimensions is not diagonal and the eigenvalues depend on 
anisotropy parameters, as well as on the parameter u. It 
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A[7,1]/£, a=a 2 , 6=3 




FIG. 34: Dependence of the critical dimension A[7, l]/f on 
anisotropy parameter ot\ = a? for different fixed point values 
of the parameter u (for notation see the caption in Fig. [5}. 
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FIG. 35: Dependence of the critical dimension A[7, 3]/£ on 
anisotropy parameter ai (ct2 — 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 



leads to the sufficient difference between anomalous di- 
mensions of the models with different time correlations 
of the velocity field. On the other hand, the fact that the 
matrix (|107f> is triangular in the isotropic case (it is also 
triangular in the case with large-scale anisotropy) is also 
important here because it allows us to assign uniquely the 



concrete critical dimension to the corresponding compos- 
ite operator even in the case with small-scale anisotropy 
and study their hierarchical structure as functions of p 
(see Ref. [27| for details). As was shown in Ref. 27] within 
the Kraichnan model, as for anomalous scaling, the lead- 
ing role is played by the operators with the most negative 
critical dimensions: for the structure functions (|25|) with 
even iV it is the operator with p = and for the struc- 
ture functions (|25|) with odd N it is the operator with 
p = 1. As we shall see, the same situation also holds in 
the general case with the finite time correlations. 



C. Anomalous scaling of the structure functions in 
one-loop approximation 

The combination of the RG representation (|77|) with 
the OPE leads to the final asymptotic expression 
for the structure functions (|25[) within the inertial range, 
namely, 



S N (r) 



JV(l-?/2) 



£ £i^'^/ i ) A[JV '" 1 + -- 



N'<N p 



(109) 

■}, 



where £ is defined in Eq. (|T§|) , p obtains all possible values 
for given N', Cm'-p are numerical coefficients which are 
functions of the parameters of the model, and dots means 
contributions by the operators others than F[N,p] (see, 
e.g., [2l|, 13 for details). 

As was already mentioned in Introduction, our aim is 
twofold. First of all, we shall find the dependence of the 
critical dimensions on the parameter u, thus we shall an- 
swer the question whether the system with finite time 
correlations of the velocity field with presence of small- 
scale anisotropy is more anomalous, i.e., whether the cor- 
responding critical dimensions are less than those of the 
Kraichnan rapid change model which was investigated in 
Ref. [13] • This question is interesting because the model 
with finite correlation time of velocity field can be con- 
sider as further step on the way to the model with veloc- 
ity field driven by the stochastic Navier-Stokes equation. 
Thus, the answer on the aforementioned question in the 
framework of the present model can also give prelimi- 
nary answer, as for possible tendencies, on the similar 
question in the framework of the scalar advection by the 
Navier-Stokes velocity field. The second aim is to inves- 
tigate whether the system with finite correlation time of 
velocity field together with the presence of small-scale 
anisotropy can lead to the more complicated structure of 
critical dimensions than it was shown in Ref. (2?| ■ There 
are two possibilities. First, it is possible that the pairs of 
complex conjugate eigenvalues of the matrix of critical di- 
mensions can exist. In this case, the oscillation behavior 
of the corresponding scaling function appears. Therefore, 
the scaling functions in Eq. (|109[) would contain terms of 
the following form 

(r/L) A « { Cl cos [A/(r/L)] + c 2 sin [A/(r/L)]}, (110) 
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FIG. 36: Dependence of the critical dimension A[7, 3]/£ on 
anisotropy parameter «2 (ai = 0) for different fixed point 
values of the parameter u (for notation see the caption in 
Fig.©. 
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FIG. 37: Dependence of the critical dimension A[7, 3]/£ on 
anisotropy parameter ot\ = oli for different fixed point values 
of the parameter u (for notation see the caption in Fig.[SJ. 



where A^j and Aj are real and imaginary part of A, and 
c\ t 2 are constants. Another, in general, possible structure 
of the matrix (| 1 07[) is related to the situation when the 
matrix of critical dimensions cannot be diagonalized and 
has only the Jordan form. Then a logarithmic correction 
would be involved to the powerlike behavior of the form 

(r/L) A [ci ln(r/i) + c 2 ] , (111) 

where A is the eigenvalue related to the Jordan cell. 

In Figs.loTI571 behaviour of the eigenvalues of the ma- 
trix of critical dimensions A[N,p] for various values of 
N = 2,3,4,5,6,7 (p = 0,2 for even values of N and 
p = 1, 3 for odd values of iV) is shown as function of the 
anisotropy parameters a\ and a 2 in three-dimensional 
case and for different fixed point values of the parame- 
ter u. The dependence of the critical dimension A[2,0] 
is not shown explicitly because it is identically equal to 
zero for all fixed point values of the parameter u. It 
can be shown either by direct calculation or by using the 
Schwinger equation (see, e.g., Ref. [13] )■ At first sight one 
can conclude that there are different behaviors of critical 
dimensions as functions of anisotropy parameters ol\ and 
OL2 and of the parameter u* for odd and even structure 
functions. Let us discuss it in detail. 

First of all we shall concentrate on the even structure 
functions (N — 2, 4, 6) and we shall discuss the behav- 
ior of the most important critical dimensions with p = 
which define the anomalous scaling of the corresponding 
structure functions. As was already mentioned in the 
case N = 2 the corresponding critical dimensions A [2, 0] 
are identically equal to zero. On the other hand, one 



can see identical qualitative behavior of the critical di- 
mensions A [4,0] and A [6,0] as functions of anisotropy 
parameters as is shown in Figs. [MllTBI and in Figs. l2"6lT2"8l 
In the case when the anisotropy parameter «2 is vanished 
the corresponding critical dimensions (as functions of pa- 
rameter ai) are the most negative in the frozen limit of 
the model (u* = 0) as is shown in Figs.Q3] and [221 On 
the other hand, in the case when the anisotropy parame- 
ter a\ is vanished (see Figs.fTSIandl!??]). as well as in the 
case when a\ = a 2 (see Figs.H6land[28)) the situation is 
opposite, namely, the most negative critical dimensions 
as functions of the corresponding anisotropy parameters 
are those that corresponds to the rapid-change model 
limit (u* — > oo). This is some kind of nonunivcrsality 
of the behavior of the critical dimensions in the plane 
of anisotropy parameters a\ — oli- Thus we still have 
some kind of hierarchical behavior in respect to u* but 
the hierarchy depends also on the values of anisotropy 
parameters. It means physically that the answer on the 
question which model is "more anomalous" can depend 
on the form of the small scale anisotropy. Besides, it is ev- 
ident that there must exist a system of curves in the plane 
a\ — «2 on which the pairs of models with different fixed 
point values of the parameter u have the same anoma- 
lous dimensions. We shall not show them explicitly here 
because we suppose that their form will strongly depend 
on the higher loop calculations which are ignored here 
(we work in one-loop approximation) but we can assume 
that the qualitative picture will be the same. Of course, 
all of the curves must cross in the point a\ — a 2 = (as 
is evident from corresponding figures for the same value 
of N) as a result of the fact that in the isotropic case the 
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critical dimensions for different values of u* are the same 
and they are given explicitly in Eq. (|108[) . 

Let us now briefly discuss the critical dimensions 
A[JV, 2] for even values of N with p = 2. Of course, 
they are not so important as critical dimensions A [AT, 0] 
but are interesting from the point of view of their non- 
trivial behavior as functions of u* as is shown in Figs. [5} 
ITI 1 TtfTQl and [291I3T1 Again one can see different behavior 
of the critical dimensions in different directions in the 
plane given by the anisotropy parameters a\ and oli but 
the most interesting feature is the fact that the corre- 
sponding couples of curves cross themselves in two points 
except for N = 2 (one of the two points is ax = a% = 0). 

As for the structure functions of odd order the situa- 
tion is slightly different. Again we start with the most 
negative critical dimensions for which p = 1. They are 
shown in Figs. lglfTUl l2~01T22l and H3J for AT = 3, AT = 5, 
and N = 7, respectively. One can see immediately that 
again in different directions in the ax — ai plane differ- 
ent models are the most anomalous (frozen limit of the 
model or rapid-change model limit), i.e., they have the 
most negative critical dimensions A [AT, 1]. But, besides, 
the situation is also different for positive and negative 
values of the anisotropy parameters. For example, in the 
case when the anisotropy parameter ai = the corre- 
sponding critical dimensions as functions of parameter 
U\ are the most negative in the frozen limit of the model 
(u* = 0) for «i > and they are the most negative in 
the rapid-change model limit of the model (u* — oo) for 
-1 < ax < as is shown in Figs.0 [2Q] and [32] On the 
other hand, in the case when the anisotropy parameter 
ax = (see Figs. [HI [22 and [33]), as well as in the case 
when ax = ai (see Figs.lTUl l2"2l and [M]) the situation is 
opposite. Thus, one can conclude that the answer on the 
question which model is more anomalous can depend on 
the form of the small scale anisotropy, i.e, on the param- 
eters of anisotropy. 

In the end, let us briefly discuss behavior of the crit- 
ical dimensions A [AT, 3] for odd values of N as shown 
in Figs.Eni [HSH and [3J][3j for N = 3, N = 5, and 
N = 7, respectively. As in the case with even values of N 
one can see different behavior of the critical dimensions in 
different directions in the plane given by the anisotropy 
parameters a\ and ai but the existence of two intersec- 
tions between couples of curves appears only for N = 7, 
i.e., it is not present in the cases with A" = 3 and A^ = 5. 

In present paper we have shown only the smallest crit- 
ical dimensions for concrete value of N, namely, p = 0, 2 
for even value of A" and p = 1, 3 for odd value of N. But 
corresponding analysis can be also done for others critical 
dimensions which correspond to higher possible values of 
p (p < N). Detail analysis shows that no exotic situa- 
tion appears in their behavior as well. Thus, we can also 
answer the second question whether the finite correlation 
time of velocity field together with small scale anisotropy 
can lead to the more complicated structure of critical di- 
mensions (oscillations or logarithmic corrections). Our 
answer is no, i.e, the matrices of critical dimensions have 



real eigenvalues at least up to Af = 7. 



VI. CONCLUSION 

Using the field theoretic RG technique and operator 
product expansion we have investigated the influence of 
uniaxial small-scale anisotropy on a passive scalar ad- 
vected by a Gaussian solenoidal velocity field with finite 
correlation time in one-loop approximation. First of all 
we have found and classified all possible scaling regimes of 
the model which are directly related to the corresponding 
IR stable fixed points of the RG equations. The "phase 
diagram" of the scaling regimes in the plane e — r\ is 
shown (see Fig. [3]) and it is found that the small-scale 
anisotropy has no influence on the stability of the scaling 
regimes (on one-loop level), i.e., we have the same five 
scaling regimes with the same regions of stability as in 
the isotropic case of the model [3(|. Two of the scaling 
regimes are related to "frozen limit" of the model, an- 
other two to the "rapid-change" model and the last one 
corresponds to general case with finite time correlations 
of the velocity field. 

Further, we have studied the influence of small-scale 
anisotropy on the anomalous scaling of the single-time 
structure functions of a passive scalar using the OPE. 
The corresponding leading composite operators with the 
smallest (the most negative) critical dimensions are stud- 
ied in detail and the critical dimensions are found as func- 
tions of the anisotropy parameters and the fixed point 
value of the parameter u which represents the ratio of 
turnover time of scalar field and velocity correlation time. 
We have shown that the corresponding anomalous dimen- 
sions, which are the same (universal) for all particular 
models with concrete value of u in the isotropic case, are 
different (nonuniversal) in the case with the presence of 
small-scale anisotropy and they are continuous functions 
of the anisotropy parameters, as well as the parameter 
u. It is shown that there is different behavior of the 
anomalous dimensions in the case of even order single- 
time structure functions than in the case of odd order 
ones, as well as there is different behavior of the anoma- 
lous dimensions in the different directions in the plane of 
the anisotropy parameters (see discussion in the end of 
the previous section for details). Thus, the answer on the 
question which special case of the general model (rapid- 
change limit or frozen limit) is more anomalous in the 
presence of anisotropy is not unique. Therefore, we are 
also not able to make definite conclusion what one can 
expect in the case of more realistic model of a passive 
scalar advection, namely, in the model of a passive scalar 
advected by the Navier-Stokes velocity field. 

It was also shown that even in the case with finite time 
correlations of the velocity field the critical dimensions of 
the corresponding composite operators have simple struc- 
ture, i.e., the matrices of the critical dimensions have real 
eigenvalues. It means that no exotic situations, namely, 
oscillations or logarithmic corrections to the critical di- 
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mensions, are present. 
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APPENDIX A 

The explicit form of the coefficients Aij with i = 1, 4 
and j = 0, 3 from Eq. (fTOTjl ) is 



Ain = 



p(p - l)[{d 2 - 5)(1 + at) + 4a 2 



d 2 - 1 



An = - 



pip - 1) 



d(d + l) 

x [d 2 + d-4-a 2 (d-7) + ai(2d 2 + d-9)], 



A i2 = P f-P lj [d + 1 + (d - l)(ai(d + 3) - 2a 2 )], 



d(d + 2) [ 



p(p - l)(a 2 - ai) 



(d+l)(d + 3) 
d(d + 2)(d + 4)' 



A 22 = -- 



1 



■{2n(ai - a 2 )(7 + 3d) 



d(d + 2) 

+n(n - 1)[12 - 60a 2 + 4ai(16 + d)] 

+np[16a 2 (d- 2) - 4ai(d - l)(d + 7) - 12(d+l)] 

+ P {p - l)[d(d + 1) - a 2 (d 2 + 2d + 9) 

+ai(3d 2 + 2d + 7)]}, 
(d + 3)(ai -a 2 ) f , , 
Ms = d(d + 2)(d + 4) i 12 ^- 1 ) 

+p(d+l)[d(p-l)-12»]}, 



2n 



1 {2(n-l)[(d 2 -l)(l + ai)-24a 2 ] 



Aso - d 2 - 
-(2p + l)(d 2 - 1)(1 + ai) - 8a 2 (d + 2 + 4p)}, 



A 3 i 



2n 



K',2 



d(d+Y) {2(n-md + l)(d + 6) 

-a 2 (d + 25) + ai(d + l)(2d + 5)] 
+ (d+l)(d + 2p(2d+l)) 
+a 2 (15 + 7d + 2p(15-d)) 
+ai(d+ l)(d(2 + 6p) - 1)}, 

In 

^^_^{2(n-l)[a 1 (18 + d(d + 13)) 

-6(a 2 - l)(d + 2)] - (d + l)[d(«! - q 2 ) 
+2(d + 2 + 3ai(d + 1) - a 2 (d + 3))p]}, 



^20 = 



-{48n(n-l)(a 2 -ai - 1) 



d 2 

-p(p - l)[16a 2 + (1 + ai)(7 + d 2 )] 
+2n[4a 2 (d + 2) + (1 + «i)(d 2 - 4d + 9) 
+2p(4a 2 + (l + ai)(d 2 -13))]}, 



.4 



21 



d(d+ 1) 

{4n(n - l)[d + 13 - 24a 2 + (d + 25)ai] 
+p(p- l)[2d 2 + d+7 + a 2 (7-d) 

+ai(3d 2 +d + 14)] 
-2n[2p(d 2 + 3d - 10) - 3d - 7 
+a 2 (2p(15-d) + 7d+15) 
+a x (2p(2d 2 + 3d - 23) + d 2 - 7d - 16)] } , 



A 
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A, 



4ra(ai - « 2 )(d + 3)((d + l)p - 6(n - 1)) 
d(d + 4) 



A iQ = -4n(n- l)(l + ai), 



\n(n - l)[2d + 4 - a 2 + 3«i(d + 1)] 



A 42 = - 



4n(n - l)[3ai(d + 2) - (a 2 - l)(d + 4)] 



.4 
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4n(n - l)(ai - a 2 )(d + 3) 
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